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Abstract 



^ ■ The dependence of the Lyapunov exponent on the closeness parameter, e, 
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in tangent bifurcation systems is investigated. We study and illustrate two 
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averaging procedures for defining Lyapunov exponents in such systems. First, 
we develop theoretical expressions for an isolated tangency channel in which 
the Lyapunov exponent is defined on single channel passes. Numerical sim- 
ulations were done to compare theory to measurement across a range of e 
values. Next, as an illustration of defining the Lyapunov exponent on many 



^ ■ channel passes, a simulation of the intermittent transition in the logistic map 

is described. The modified theory for the channels is explained and a simple 
model for the gate entrance rates is constructed. An important correction 
due to the discrete nature of the iterative flow is identified and incorporated 
in an improved model. Realistic fits to the data were made for the Lyapunov 
exponents from the logistic gate and from the full simulation. A number of 
additional corrections which could improve the treatment of the gates are 
identified and briefly discussed. 



1. Introduction and Background 



Chaos is the study of dynamical systems that have a sensitive dependence on initial 
conditions. Much attention has been paid to the two main routes to chaos: pitchfork 
bifurcation and tangent bifurcation. If we consider the general difference equation mapping, 

x n+1 = F(x n ), (1) 

then tangent bifurcation, also called type I intermitency [Pomeau & Manneville, 1980], 
occurs when a tangency develops in iterates of F(x n ) across the x n = x n+ i reflection line. 
(Pitchfork bifurcations occur when iterates of F(x n ) possess perpendicular crossings of this 
line.) Just before the tangency occurs (characterized by the closeness parameter, e, being 
small), the map is almost tangent to the reflection line and a long channel is formed. When 
the iterations enter, a long laminar-like flow is established, with nearly periodic behavior. 
Once the iterations leave the channel, they behave chaotically, then re-enter the channel. 
The result is a long region of laminar flow that is intermittently interrupted by chaotic 
intervals. This occurs when e is near zero and tangency is about to occur, hence the two 
names: intermittent chaos and tangent bifurcation. Experimentally, type I intermittency 
has been observed in turbulent fluids [Berge et. al, 1980], nonlinear oscillators [Jeffries & 
Perez, 1982], chemical reactions [Pomeau et. al., 1981], and Josephson junctions [Weh & Kao, 
1983]. An excellent introduction to the intermittency route to chaos is given in Schuster 
[1995]. 

In the pioneering studies [Manneville & Pomeau, 1979] and [Pomeau & Manneville, 
1980], it was found that the number of iterations followed an e' 1 ^ 2 dependence and that the 
Lyapunov exponent varied as e 1 / 2 for a logistic mapping [z — 2). In the work by [Hirsch et. 
al., 1982], an expression for the number of iterations spent inside the channel was developed. 
The equation for the third iterate, i.e. F(F(F(x))) or F^(x) where F(x) = Rx(l — x), was 
expanded in a Taylor series about one of the tangency points for R c — 1 + y/8. In the case 
of the logistic map, we get 
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F( 3 \x) =x c + (x- x c ) + a c (x - x c ) 2 + b c (R c - R), (2) 

where x c is one of the three contact points. After a transformation that centers and rescales 
the system around x c (y n = Xn ^ Xc ), the recursion relation can be put into the form 

y n+1 = ay 2 n + y n + e, (3) 

where e = R c — R > and a = a c b c . The more general case can be studied as a first, second, 
or any iterate instead of just the third iterate, as long as a tangency develops. To derive 
an analytic description of the trajectory, [Hirsch et. al, 1982] switched from a difference 
equation to a differential equation. Thus, they considered 

|=^ + e . (4) 

This approximation is justified as long as the number of iterations in the channel is large 
enough or, alternately, that the step size between iterations is small compared to the channel 
length. This is an easy differential equation to solve. One obtains 
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u yin is the entrance to the tangency channel and u y out " is the exit value and one has that 

- Vout < Vin < Vout- (6) 

[Hirsch et. al, 1982] observed that the entrance points for the logistic map, y in (R c > R), 
had a probability distribution that was roughly uniform. Given this distribution, the average 
number of iterations to travel the length of the channel is given as 

<n>= / n(y in )dy in = —= tan [y ut\ - ■ (7) 



[Hirsch et. a/., 1982] also derived a form for the average number of iterations for an arbitrary 
universality class. The universality class, z, is given by the lowest non vanishing power of 
(x — x c ) in the expansion around the tangency point. For tangency to develop, z must always 
be an even number: 



y n+1 = ay z n + y n + e. 



(8) 



This leads to the differential equation, 
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and to the number of iterations, 
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The average number of iterations is given by 




dy 
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(11) 




when the entrance distribution is again uniform. The numerical simulations in [Hirsch et. 
al, 1982] agreed with predicted values quite well. 

There are two manners in which Lyanunov exponents may be defined in a simulation with 
many trajectories. One may define a procedure which measures the Lyanunov exponent on a 
given trajectory, for example a single channel pass, and then averages over these trajectories. 
Another possibility is to measure the exponent across many trajectories or channel passes, 
using a binning procedure to define variances. We will use both procedures here to illustrate 
the theory. The first procedure will be termed a single pass measurement, the second a many 
pass measurement. We will develop the theory for the first procedure in the next Section, 
which will then be illustrated in Section 3 by a simulation in an isolated tangency channel for 
general z. As an illustration of a many pass measurement, a simulation of the intermittent 
transition in the logistic map will be described in Section 4. The modified theory will be 
motivated and a simple phenomenological model of the data will then be given in Section 5. 
In Section 6 an improved expression for the inverse number density, ^ due to the discrete 
nature of the iterative flow will be developed. This will improve the comparison of the model 
to measurement. Finally, we will summarize our findings and make suggestions for further 
improvements in the model in the final Section. 
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2. Tangency Channel Theory 



Our analysis of the system described by Eqs.(|8|) and @ is built on the work of both 
[Pomeau and Manneville, 1980] and [Hirsch et. al, 1982]. In contrast to the situation 
for the average number, < n >, little work has been done to develop expressions for the 
Lyanupov exponents for the tangency channel in intermittent systems. We are interested in 
understanding the e dependence of the Lyapunov exponent for z — 2, finding the constant 
of proportionality, and generalizing the results for an arbitrary universality class, z. 

The Lyapunov exponent is a measurement which characterizes the sensitive dependence 
on initial conditions of chaotic systems. It is defined as the coefficient of the average expo- 
nential growth per unit time between initial and final states of a system, which in this case 
we will take to be a single pass through a tangency channel. It is given in the case of the 
one-dimensional mappings considered here by [Scheck, 1994] 

A S ]imi£.n« (12) 

This gives us our starting point for deriving a theory for the Lyapunov exponent for a system 
with an arbitrary universality class. In that case, the function F(yi) from (|8|) is 

F(Vi) = ay* + Vi + e, (13) 

so 

Since we are interested in the Lyapunov exponent for the tangency channel, there are only 
a finite number of steps during which the iterations are confined to the channel and the 
appropriate value for n for this trajectory is the total number of iterations in the channel, 
n{yi n ). With this in mind, the Lyapunov exponent is modeled by 

^ Vout 

KVin) = — r — r / dn In |1 + azy z ~ l |, (15) 



where we have replaced the discrete sum by an integral over n-space. Again, this step is 

justified as long as the number of iterations is large enough so that the values of the natural 
logs of the slope are almost continuous. 
From Eq.(|§D, we have 

dn = — , (16) 

ay z + e v ' 

so that 

^ = S£j J dV ov + e ■ (17) 

Vin 

This gives the Lyapunov exponent for the system starting at yi n and ending at y out . 

In the logistic map or any other system, the entrance into the tube is random. Since 
the starting points are randomly distributed, it is more useful to derive a formula for the 
average Lyapunov exponent per pass. Using the above formula for \(yi n ), we can calculate 
the average value of the Lyapunov exponent and obtain 

Vout 

< A >= J dy in X(y in )P(y in ), (18) 

—yout 

where the probability function P(yi n ) satisfies 

yout 

J dy m P(y m ) = l. (19) 

— yout 

For the present, let us consider the special case of a uniform distribution, 

PiVin) = (20) 

Using this probability distribution, we obtain 

yout . . 

<A>=^ / dy in ^{, (21) 
2y<mt J n(y in ) 
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where 
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and where n(yi n ) is given by Eq.flTOD above. 

One approximation and a change of variables are necessary to make this formula more 
usable. One important step is to define the value for y out as 

Vout = s{ -, (23) 
V a 

where "s" is a positive scale factor that can be independently set in order to model a given 
system. Clearly, s can not be arbitrarily large. A natural requirement is that the derivative 
in ([y]) be positive, making the possible "throat" of the channel end at the point where 
= 0. This gives that 

Smax = z^a^^e' 1 -, (24) 

is the maximum value of s for given z, e and a. With the value of y ou t from (|23"D , the integral 
I{yin) can be simplified with a change of variables. Let 

y' = —■ (25) 

Vout 

Therefore 
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I{y in ) = sa-*e-^ J dy'^ — — (26) 

where the absolute value in the natural log is no longer necessary. The Taylor series expan- 
sion for natural log is (\x\ < 1) 

2 3 4 
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Hl + x)=x- _ + _-_ + ... . (27) 
Using this approximation we have simply 

as long as 

s << s max . (29) 



Hv*)***' I dy-S^ = ^(^ 1 r)^ (28) 



Our simplified formula for the average Lyapunov exponent is now 



< A >= U 1 '^ 1 ' 7 In f £±i_) ■ (30) 



—yout J y 2 +l 

Vin yf 

We now make the same scale change in the ?/j„ integral as in the y integral in (p5*|): 

^ (31) 

Vout 

Therefore 

< A >= Vv-V' / ln f + * ) . (32) 
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As one can see, for constant s the average Lyapunov exponent varies as e 1-1 ^ with a constant 
of proportionality determined by the parameters a and s. In the case where z = 2 this gives 

<\>=-^a~e ( T —^- t- — — ln ( S [ + 1 ). (33) 

2 V J tan' 1 (s)- tan' 1 (sy) \s z y z + l) y ' 

For a general probability distribution, we would have instead 

< A >= a^e 1 ^ j dy^~ ln ( dy, (34) 

sy 

where 

J 1 i dyP(y) = l. (35) 

In the case of constant scale factor s, we therefore see that the single pass tangency channel 
Lyapunov exponent behaves like e 1_1//z . 

3. Tangency Channel Simulation 

The integral in Eq. (|33]) was calculated using numerical methods and compared against 
numerical simulations of the logistic map (z — 2). In all cases we used a simulation consisting 
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of 10,000 Monte Carlo runs for each data point. As can be seen in Figs. 1-3, for values of 
s = 0.1, 1.0 and 10, the theoretical values agree with the simulation values for a uniform 
probability distribution for small enough e. At low e, the agreement is excellent, with 
the theoretical value straddled by the upper and lower error values of the simulation. In 
the s = 0.1 simulation, the assumption that the discrete sum can be approximated by a 
continuous integral breaks down at large enough e due to the very small number of iterations 
in the channel. For the s = 10 simulation, the natural log approximation, Eq . fl27D , starts 
to break down and is the main cause of the divergence between theory and simulation. 
The least divergence between theory and simulation occurs when s ~ 1 . The calculations 
become more time consuming at larger s due to the increased number of iterations necessary 
to pass through the channel. 

We examined the more general expression, Eq.(|4]), when z = 2 for other probability 
distributions including a Gaussian and a \y\ distribution of normal deviates. Although these 
results are not illustrated here, the theoretical and simulation results were again in excellent 
agreement for small enough e. We also examined the e dependence for higher universality 
classes. However, due to the large amount of computer time it takes to do such simulations, 
we have data only for one additional z value. For a universality class of z = 4, the Lyapunov 
dependence should be e 3//4 , which is clearly confirmed in Fig. 4. 

The Section 2 expressions for < A > are displayed in a form appropriate for a simulation 
at a fixed value of the scale, s. However, the approximation employed, Eq. fl29|) , is simply 
the condition that the channel length be much smaller than the total gate size (determined 
by the tangency point and the point at ^ = 0). Thus, the expression Eq.(|34]) holds also for 
a fixed tangency channel fraction, /, 



./ 



Vout 



(36) 



az) 1 ~ z 



as long as 



/«!■ 



(37) 
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The quantity (az) is just the total model gate size. For z = 2 the relationship between / 
and s is simply 

/ = 2s^a~e. (38) 

When the change from s to / is made in Eq . fl34|) , the result is no longer proportional 
to e 1 / 2 but in fact goes to a constant at small values of e. Mathematically, this is due to 
the fact that the denominator, proportional to the number of iterates in the channel for a 
starting position y, falls off like e 1//2 when y ~ 1. Physically, the Lyapunov exponent is being 
dominated by the small number of iterates associated with entrances on the far side of the 
narrow channel. Fig. 5 shows the predicted and measured values of the Lyapunov exponent, 
< A > for the case / = 0.1, i.e., the channel length is one-tenth the size of the gate, when the 
entrance probability is again uniform. As expected, and unlike the cases presented above 
at fixed s, the value of the Lyapunov exponent becomes constant at small e, the theoretical 
value remaining about 10% larger than simulation. This amount of deviation is what we 
would expect since the term kept in the expansion of the natural log in Eq.([27|) is of order 
/ and the first neglected term is of order f 2 . Measurements at fixed / will be important for 
the simulations done in the full logistic map to be considered in the next Section. 

4. Logistic Map Simulation 

As was pointed out earlier, there are two manners in which Lyapunov exponents can be 
defined in channel simulations, which we called the single pass and many pass definitions. 
We have already described and illustrated a single pass simulation. As an example of a many 
pass situation, we consider a simulation of the tangency gates in the logistic map. In doing 
so, we find it convenient to first describe the details and results of the numerical simulation. 
The theory and the model used to fit to the data will then be developed together in Section 
5. 

As outlined in the Introduction, the equation for the third iterate of the logistic map for 
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small e may be expanded in a Taylor series about each of the tangency points. The positions 
of the tangency points are given to high precision in Table I. (Knowledge of any one gives 
the others through the basic recurrence relation.) Expanding to third order in (x — x c ), one 
has for the third iterate, 

F {3 \x) =x c + (x- x c ) + a c (x - x c f + c c (x - x c f + b c (R c - R). (39) 

The values of a c , b c and c c are given in Table I. Again introducing y n = Xn ~ Xc ; the difference 
equation for all three of the gates takes the form (note that c c b 2 c = —196 to high accuracy 
for all gates), 

2 

y n+ i -y n = e + ay 2 n - ^y^, (40) 

where the constant "a "takes on the value 

a = 69.29646455628 .... (41) 

An interesting aspect of the simulation is the exclusion of certain x-values from the 
logistic map at finite e. We have labeled the tangency gates in increasing order of their x c 
values. Referring to the third iterate map, shown for e = in Fig. 6, it is clear by drawing a 
horizontal line that gate 1 is reachable only under steady-state conditions from points close 
to point C in the Figure. Likewise, gate 3 is reachable only from previous iterates starting 
on or near points A or B in the figure. (Since x = 0.5 is a symmetry point in the mapping, 
the values of F^ 3 \x) at points A and B are the same.) Note in this context that the laminate 
flow through gates 1 and 2 is from smaller x-values to larger ones, whereas the flow through 
gate 3 is in the opposite direction. Iterates entering gate 1 from the point C will actually 
enter at the value xl = F^ 3 \xc)] points between x = and x = xl will never be reached. 
This is after a possible transient of a single iteration. Likewise, iterates entering gate 3 from 
points A and B will enter at the value xr = F ( - 3 \x a,b)] points between x = xr and x — 1 
are never reached, again after a possible 1-iteration transient. By drawing a horizontal line, 
gate 2 is seen to be reachable from 4 separate x-regions (excluding the points to the left of 
Xl and to the right of x R ). 
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The measurements leading to values and variances of the Lyapunov exponents were taken 
from a single trajectory of 1.6 million iterations of the third-iterate logistic map at each e 1 / 2 
value following an initial "heating" to remove the transient x- values. Monte Carlo error bars 
for the Lyapunov exponents were then measured by breaking the single trajectory into 100 
bins. Runs were made both at fixed gate fraction, /, as well as at fixed scale factor, s. 

Fig. 7 shows a log 10 plot of the contribution of each of the three gates when the fraction 
of the gate being measured is / = 0.1, the same used in Fig. 5. Note that / = 0.1 indicates 
the fraction of the model gate size, characterized by Eq. (H), not the actual gate size. The 
ratio of the model to actual gate length is about 1.025. The gate Lyapunov exponents, 

< A > g , are normalized to the number of third-iterate hits within each gate, which at small 
e just approaches 1/3 of the total iterations. We notice that even at the larger e values the 
contribution from different gates is the same within errors. 

Fig. 8 shows the contribution of each of the three gates to the Lyapunov exponent, 

< A > g , in a simulation with a fixed value of the scale factor, s = 1.0. Eq.(|7|) implies that 
such a simulation will have a mixture of half iterations inside the gates and half outside, as 
indeed is observed. The error bars here are larger in a relative sense than in Fig. 7 because 
the gate size is shrinking like e 1//2 (see Eq.(^3j)), leading to smaller statistics. Also unlike 
Fig. 7 there are significantly different contributions from the three gates at larger e values, 
the middle gate having an enhanced < A > g . It is only at values of e below and including 
e 1 / 2 = 0.512 x 10~ 3 that a distinction between the gates can no longer be seen. However, 
this may simply be the result of the larger statistical fluctuations present at smaller e. 

Note that the fixed / data in Fig. 7 shows an approximate e 1//2 behavior, while the fixed 
s simulation in Fig. 8 behaves approximately as e at small e values. These behaviors are 
in contrast to the simulations in Section 3 where the fixed / data went to a constant at 
small e and the fixed s data behaved like e 1 / 2 . The e behavior of the logistic map Lyapunov 
exponents will be commented on further in the next Section. 

In order to model the Lyapunov exponents for the gates, it is necessary to have an 
understanding of the entrance probability for the gates as a function of position in the gate. 

11 



There are two ways in which a laminate flow may begin in the tangency channel. Primarily, 
entrance into the gate occurs as a continuous flow from just outside the gate. Alternatively, 
the flow can begin in a discontinuous fashion from a disjoint region of the map. These two 
entrance routes will be termed the continuous and discontinuous types, respectively. Fig. 9 
presents a measurement of the binned discontinuous entrance rate, n^, for the first gate of 
the / = 0.1 simulation at e 1 / 2 = 0.128 x 10 -3 . The data is divided into 19 bins with bin 
size of Axbin = 0.5624 x 10~ 4 , which is just the model gate width divided by 100. The first 
bin (which would have extended from -10 to -9 in the figure units) contains both continuous 
and discontinuous entrances. Since we have not attempted to separate the discontinuous 
from the continuous entrances in this bin and because the total rate is off scale, this first 
bin is not shown. The entrance rate seems to be fairly uniform in this Figure; gate 2 and 
(the mirror image of) gate 3 look very similar. This approximate entrance uniformity for 
small / will be useful in setting up a simple model of the gate contribution to the Lyapunov 
exponent, which will be described in the next Section. 



We will now develop theoretical expressions for the e-dependence of the Lyapunov expo- 
nents for the logistic map. For this purpose we need to develop expressions for < n > and 
< A > in a many pass simulation, as opposed to the single pass considerations in Sections 
2 and 3. In a many pass simulation, the number of iterations in the gate will be weighted 
by the entrance rate rather than probability. Therefore, Eq.(^), generalized to an arbitrary 
probability distribution, is replaced with 



5. Simple Model 




(42) 



where 




(43) 
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We also are using the dimensionless variable y introduced in Eq. (|3~TD (y,„ = yy ou t)- The 
functional form of dN d ^ has yet to be specified. Note that nd in Fig. 9 is given by nd = 
^^Ay bin where in this case Ay bin = 0.1. 

The form for the Lyapunov exponent, Eqs. (|2l|) and (|22|), must also be modified. It is 
again the rate rather than the probability which is relevant. In addition, for a many pass 
simulation the Lyapunov integrand in Eq. (|2TD must be weighted by the ratio of the number 
of iterations for each passage through the channel, n(yi n ), to the average total number of 
iterations, < n > g , resulting in 

< A ^2^/>^'<^ < 44 > 



where 



2 ri , ln(l + fy') 



I(Vin) = 7 / dy' /2 4oe ■ ( 45 ) 

J J y y t jt 
Making the same approximation as in Eq.([27]) above, this simplifies to 



1 f 1 ,~dN(y)^ ( 1 



2 



< x >^2^L di ^{WT^^ (46) 

We use / instead of s in these formulas since most of the simulations in the following will 
use fixed /. 

As we have seen in the last Section, the Lyapunov exponents from the three gates are 
apparently indistinguishable at small enough e. A very useful simplification therefore is 
to ignore the distinction between the gates and model the Lyapunov exponent as if there 
were only two regions, the gate (or periodic) region and the outside (or chaotic) region. In 
addition, Fig. 9 suggests that a reasonable model of the channel region is to assume that 
the discontinuous entrance rate is uniform. This last simplification is only possible for small 
enough /, the entrance rate in the complete tangency channels being far from uniform. 
These assumptions will allow us to construct a very simple model of the e dependence of 
the Lyapunov exponents. The choice of / = 0.1 to separate the two regions is arbitrary. 
A smaller value would result in an even more uniform entrance rate than Fig. 9; however, 
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one would also loose statistics because of the smaller gate size. As we will see, the gate 
fraction / will be used to formally separate the outside and inside gate Lyapunov exponent 
behaviors. 

A new aspect of modeling the actual gates of the logistic map is the fact that the entrance 
to the gates can be from flow further up the channel or from a completely disjoint part of 
the map. These possibilities were termed the continuous and discontinuous routes in the 
last Section. One may show that the continuous entrances occur within a scaled distance 
of ~ f /2 + 2ae/f from the lower limit (—1) of the integrals in Eqs. fl42|) and (f44|). In Fig. 9 
this contribution would extend about halfway through the first (deleted) bin. Since these 
entrances always occur in a narrow range of the integrations in these equations for the range 
of e considered, it is reasonable to model this contribution by a Dirac delta function located 
at the lower limit, y = —1. In addition, we saw in Fig. 9 that the discontinuous part of the 
entrance rate was approximately uniform. Thus we will model the entrance rate with 

dN(y) OAr , / , dN d 

= 2N c 5{y + 1) + -jJL, (47) 

where N c and are constants in y. This gives a simple model of the contributions to 
< n > g and < A > g from the logistic gates: 

< n > = _L tan-^-^) (2N C + , (48) 

yjae 2y/ae \ ay J 



< A > 9 = ^ t dyln ] . (49) 

9 2 < n > g dy J-i \y 2 + ^£ ] 



Notice that N c drops out of the expression for < A > g . In fitting the data, we also need the 
expression for the small e limit of Eq.(f48l), which we will call n^: 

7T / dNA , . 

As pointed out previously, / represents a separation parameter for the gate and outside 
regions. Inside the gates one expects from the previous results that the number of iterations 
associated with a given traverse of the gate will increase like e -1 / 2 at small e, independent of 
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the form of the entrance probability. Thus in the many pass simulations of the full logistic 
map, one expects any given fraction / of the gates at small e to eventually contain essentially 
all iterates. This means from Eq.(|50|) that the quantity 2N C + must scale like e 1 / 2 for 
ut to become constant. This behavior will be assumed for these quantities individually. 
In addition, we assume that for small enough gate fraction the continuous entrance rate is 
uniform. Given these assumptions, the parameters N c and can be parameterized as, 



N c = v^(l - f)n T K c , (51) 
dN d 



= V^fn T K d , (52) 

where K c and K d are constants. 

We can now better understand the e dependencies in the Lyapunov exponents in the 
logistic map simulation seen in the last Section. For a many pass simulation we expect the 
model Lyapunov exponent from Eqs. (f49|) and ( 52 ) for fixed / to behave like e 1//2 at small e. 



This is because < n > g saturates to the value tlx while the rate in Eq.(^) goes like e 1 / 2 . 
In contrast, for fixed scale s the exponent now acquires an extra e 1//2 factor from the gate 
fraction / in ( [52] ) and is expected to decrease like e, as was seen in the actual simulations. 
This is just a result of the shrinking gate size of the fixed s simulation. 

As a result of the number flow into the gate regions as e decreases, the outside region 
becomes sparsely visited, but with the same local density of iterates. (The shape of the 
outside region changes very little for small e.) This implies the Lyapunov exponent measured 
in the outside region will go to a constant for small e. The flows being described are all due 
to the result Eq. (pT0|) for n(y in ) and can be thought of as applications of the renormalization 
flow arguments, without external noise, in [J. E. Hirsch et. ai, 1982] and [ B. Hu and J. 
Rudnick, 1982]. 

Although the emphasis here is on the gate Lyapunov exponents, one can now make a 
rough model of the complete logistic map exponent. Letting < A > g represent the Lyapunov 
exponent expression for the fixed / gate (laminar region) from Eq.(^9[) and < A > Q be the 
outside (chaotic region) contribution, the expression for the Lyapunov exponent for the 
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complete logistic map in this model is 

< A > 3rd = — < \ > < A > g , (53) 
n T n T 

where n a + n g = rix- That is, < A > 3rd is just assumed to be a sum of the exponents 
< A > g and < A > Q weighted by the relative number of iterations spent in the two regions. 
A more fundamental description of the logistic map would independently calculate < A > Q . 
However, from the previous arguments we expect < A > to simply be a constant at small 
e. It will be evaluated from a fit to the data. 

Numerically, the gate contribution in Eq . (|53|) is only about 1% of the total through 
most of the e range considered for / = 0.1. Both terms in Eq . (|53"D go like e 1 / 2 , but in 
different ways. The outside term behaves like e 1//2 because the outside number n has this 
dependence; the gate term also behaves this way because < A > g itself goes like e 1 / 2 at fixed 
/ as explained above. Note that all Lyapunov exponents in Eq.([)3|) are normalized to the 
number of third-iterate steps in the simulation. This is symbolized by writing < A > 3rd for 
the complete logistic map Lyapunov exponent. We must remember to divide this value by 
three to calculate the usual single-iterate value: 

< A > lst = - < A > 3rd . (54) 



Fig. 10 presents the results of fitting Eq.fl49|) to the data for the / = 0.1 gate exponent 
and Fig. 11 gives the measured and model < A > lst values for the complete logistic map. 
There are three parameters needed to do these fits: K c , K d and < A > Q . To evalute these 
constants we simply fit the values of these expressions to the measured values of < A > ls * 
and < A > g at a value of e 1 / 2 = 0.128 x 10~ 3 , near the middle of the exponential range in 
these figures, as well as the maximum number of gate iterations, ut- Since we are averaging 
over the properties of all three gates, rir — § x 1.6 x 10 6 for the simulation in Section 4. 
This gives K c = 0.358 x 10" 2 , K d = 0.739 x 10~ 2 and < A > Q = 0.962. 

Of course since Figs. 10 and 11 were used to fit the model parameters, we need an 
independent test of how well the model truly represents the data. For this purpose we also 
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present Fig. 12 and Table II. Fig. 12 compares the model results against the data for a 
s = 1.0 simulation. As explained above, this data falls like e. The theoretical results are 
satisfactory although seem somewhat high compared to measurement. In addition in Table 
II we give the fit results for the rates N c and compared to measurement. The measured 
value for ^¥ is actually just an average over all three gates of binned data similar to Fig. 9, 
and the value for N c is the average value in the three entrance bins minus the average of 
the entrances in the other bins. The results for iV c are good but the fit values of are 
approximately a factor of 2 larger than measurement. As we will see in the next Section, 
this is largely the result of an inaccurate characterization of the inverse number density, j*-. 



6. Improved Model 

Although the expression for the inverse number density, Eq.([|), is symmetric (even) in 
y for z — 2, there are at least three sources of asymmetry in the actual gate. First, the 
discontinuous entrance rate raises the value of the hit density on the exit sides of all 
three gates. Second, the term proportional to (x — x c ) 3 in Eq.(^) shows that there is a 
small intrinsic asymmetry in the shape of the gates themselves, raising the hit density in 
the same sense. Most interestingly, there is also a contribution due to the finite step size of 
the laminar flow through the gates which contributes even for a perfectly symmetric gate. 
This will now be described. 

Remembering that a finite difference gave rise to the left hand side of Eq.(|]), a more 

accurate differential characterization of the iterative flow is 

dy ld 2 y 2 , , 

d^ + 2d^ = ay+e - (55) 

We will use the method of successive approximants to evaluate the second derivative term. 

To zeroth order, 

_| = ay + e. (56) 
Thus the lowest order result for the second derivative is just 
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2ay(ay 2 + e). 



(57) 



dn 2 



o — 



Inserting this back into the starting point, Eq.(^), we now have the improved result 



dy_ 

dn 



i = ay 2 + e(l - ay) - a 2 y 3 . 



(58) 



With this improvement, a better formula for the gate Lyapunov exponent is given by 



Note that the intrinsic contribution to the inverse hit rate from Eq. ([40"D is of the same sign 
but 49/2 times smaller than the term from Eq . fl58|) and so is neglected. (The intrinsic term 
would also slightly alter the numerator; see Eq.fljjD.) 

Notice that the continuous entrance term, N c , now does contribute to the expression for 
< A > g since the inverse number density is no longer an even function. Unfortunately, the 
innner integral can no longer be done analytically and a double integral survives. Because 
the emphasis here is on modeling the Lyapunov exponent and because of the difficulty 
of performing the numerical integration leading to tit at small values of e, we have not 
attempted to make the same correction in the expression for < n > g . Thus, we continue to 
use the expression Eq.(|5|) above for < n > g in the gate region. 

When the same sort of fit is made to the simulation data as in Section 5, there is sur- 
prisingly little change in the functional forms in Figs. 10, 11 and 12, although the inverse 
number rates of the two models are considerably different and the N c term is now contribut- 
ing about 50% of the total. The new fit gives K c ~K d = 0.379 x 10~ 2 . (The value of < A > Q 
does not change from the previous fit.) The major improvement occurs in the value of the 
discontinuous entrance density, ^jf (see the "Improved model" columns of Table II), which is 
now within 5% of the measured value at e 1 / 2 = 0.128 x 10~ 3 , where the fit is actually made. 
However, the value for N c has increased and is now approximately 6% large compared to the 
measured value. (Note that 2N C + ^4 in Table II is required to have the same value in the 




(59) 
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two models because the form for < n > g is unchanged.) This problem should be cured when 
the more accurate result for the number density implied by Eq. (|58]) is used in the expression 
for n(y in ) in Eq.(|H]). 

7. Summary and Conclusions 

Tangent bifurcation or intermittent chaos is a common occurrence in systems that exhibit 
chaotic behavior. In these systems the intermittent behavior can be modeled by differential 
and difference equations of some universality class. We found that the Lyapunov exponent 
for isolated gates at single channel pass can be modeled given the universality class, the 
parameters of the difference equation, the scale factor s or fraction / of the gate size, and 
the closeness factor e. Our main theoretical result for these systems, subject to the restriction 
of a sufficient number of steps in the channel and the small gate approximation in Eqs.(^) 
or (|37|) , is that the average Lyapunov exponent is given by Eqs.(|4]) and (|35|). Single pass 
numerical simulations were consistent with these expressions. At fixed scale factor s these 
results gave a Lyapunov exponent proportional to e 1_1//z for a tangency channel with general 
universality class, z. We also showed that a simulation at fixed gate fraction / gave a result 
which instead became constant at small values of e for z = 2 due to a small number of 
entrances on the far side of the narrow channel. 

Simulations were also performed on the full logistic map near the intermittent transition 
at R = 1 + y/8. Modified expressions for the gate number, Eq. fl42|) , and Lyapunov exponent, 
Eq. (fj6|) , for a many channel pass simulation were motivated. A new aspect encountered 
in the description of the actual tangency channels was the existence of a continuous flow 
contribution into the tangency gate. Two phenomenological models of the channel were 
constructed and examined. A very simple model was considered which was able to give a 
fairly realistic characterization of the various Lyapunov exponents and the continuous, N c , 
and discontinuous, entrance parameters. We also derived a first-order correction to the 
inverse hit rate due to the discrete nature of the iterative flow, which mainly improved the 
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comparison with the measured discontinuous entrance parameter. 

Besides applying the finite discretization correction to the number density expression, 
there is considerable room for improving the present model of the logistic gates. For exam- 
ple, the modeling of the continuous contribution to the gates as a delta function is clearly 
oversimplified; no attempt has made no attempt to resolve the shape of the continuous en- 
trance rate in the binning procedure of Fig. 9. In addition, the approximation used for the 
logarithm, Eq . fl27|) , can be removed at the cost of a more complicated numerical evaluation. 
Finally, further discretization corrections to both < A > g and < n > g should result in an 
improved characterization of the inverse number density, leading to better comparison with 
the measured rates and functional behaviors at larger e values. 
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TABLES 



TABLE I. Values of the tangency points, x c , and the constants a c , b c and c c in Eq.(39). 





x c 


a c 


b c 


c c 


gate 


1 


0.1599288184463... 


88.91012989368... 


0.7793989800616 . . . 


-322.6535182739... 


gate 


2 


0.5143552770620... 


34.14530797001... 


2.029457886780 . . . 


-47.58783903539 . . . 


gate 


3 


0.9563178419736 . . . 


-310.6483669763 . . . 


-0.2230704292148... 


-3938.873792041 . . . 



TABLE II. Results from two models for the continuum number contribution, N c , and the 
discontinuous number density, % when / = 0.1. 





Simulation 


Simple 


model 




Improved model 




N c 


dNd 
dy 


N c 


dNd 
dy 


N c 


dNd 

dy 


0.8 x 10~ 5 


0.1021(6) x 10 2 


0.112(5) x 10 1 


0.101 


X 


10 2 


0.232 x 


10 1 


0.107 x 


10 1 


0.119 x 


10 1 


0.16 x 10" 4 


0.2028(8) x 10 2 


0.224(6) x 10 1 


0.203 


X 


10 2 


0.464 x 


10 1 


0.214 x 


10 2 


0.237 x 


10 1 


0.32 x 10~ 4 


0.403(1) x 10 2 


0.466(8) x 10 1 


0.406 


X 


10 2 


0.927 x 


10 1 


0.429 x 


10 2 


0.474 x 


10 1 


0.64 x 10" 4 


0.804(2) x 10 2 


0.93(2) x 10 1 


0.812 


X 


10 2 


0.185 x 


10 2 


0.857 x 


10 2 


0.949 x 


10 1 


0.128 x 10~ 3 


0.1607(2) x 10 3 


0.181(2) x 10 2 


0.162 


X 


10 3 


0.371 x 


10 2 


0.171 x 


10 3 


0.190 x 


10 2 


0.256 x 10~ 3 


0.3181(3) x 10 3 


0.368(2) x 10 2 


0.325 


X 


10 3 


0.742 x 


10 2 


0.343 x 


10 3 


0.380 x 


10 2 


0.512 x 10~ 3 


0.6274(4) x 10 3 


0.718(4) x 10 2 


0.649 


X 


10 3 


0.148 x 


10 3 


0.686 x 


10 3 


0.759 x 


10 2 


0.1024 x 10~ 2 


0.1216(1) x 10 4 


0.1411(5) x 10 3 


0.130 


X 


10 4 


0.297 x 


10 3 


0.137 x 


10 4 


0.152 x 


10 3 


0.2048 x 10~ 2 


0.2297(1) x 10 4 


0.2623(7) x 10 3 


0.260 


X 


10 4 


0.594 x 


10 3 


0.274 x 


10 4 


0.304 x 


10 3 


0.4096 x 10~ 2 


0.4135(2) x 10 4 


0.4769(9) x 10 3 


0.519 


X 


10 4 


0.119 x 


10 4 


0.548 x 


10 4 


0.607 x 


10 3 
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Figure Captions 



1. Simulation of the system described by Eq.(^) compared to the Lyapunov exponent 
given by Eq. (|33|) (z = 2; uniform entrance probability). We plot log 10 < A > against 
log 10 e 1//2 ; the prediction ([!3|) is given by the dotted line. The Monte Carlo error bars 
on the calculation are extremely small and are given by the data point bars. We are 
using a = 34 (the same as in [Hirsch et. al, 1982]) and s = 0.1, with e 1//2 ranging in 
value from 0.8 x 10~ 5 upwards by factors of 2. 

2. The same as Fig. 1 but for s=l. 

3. The same as in Fig. 1 except for s = 10. Note that the largest e 1//2 value, present in 
Figs. 1 and 2, violates the bound s < s max of the text and has been excluded. 

4. The case of z — 4, s = 0.1 a = 34 and uniform entrance probability. 

5. Contribution to the Lyapunov exponent, as a function of log 10 e 1//2 , from a simulation 
involving 10,000 gate entrances when /, the gate fraction, is set to 0.1, the entrance 
probability is uniform and a = 34. The theoretical result from Eq.(|3~3"|) is given by the 
dotted line. 

6. The third iterate of the logistic equation, F^ 3 '(x), as a function of x when e = 0. The 
three points where the third iterate makes tangential contact with the 45-degree line 
are the tangency points. The points A, B and C in the map are discussed in the text. 

7. Contributions of the three tangency gates to the logistic map Lyapunov exponent, 

< A > g for / = 0.1. Gate 1 data is given by the circles, gate 2 by the squares and gate 
3 by the triangles. Note that the ordinate values of the gate 1 and 3 data points have 
been shifted to the left and right, respectively, for clarity of presentation. 

8. Contributions of the three tangency gates to the logistic map Lyapunov exponent, 

< A > g for s = 1.0. The meaning of the symbols is the same as in Fig. 7. 

23 



9. The number of discontinuous entrances, n d) for gate 1 when / = 0.1 and e 1 / 2 = 
0.128 x 10 -3 as a function of bin number centered about the first tangency point, x c . 
Each data point above corresponds to entrances in a bin size of Axu n = 0.5624 x 10~ 4 . 
(The bin extending from -10 to -9 is not shown; see the text.) 

10. Comparison of the model results for the Lyapunov exponent for the / = 0.1 simulation 
with the averaged data from Fig. 7. The model values, given by a dotted line (simple 
model) and a solid line (improved model), are indistinguishable. The data point at 
log 10 (0.128 x 10~ 3 ) = -3.893 ... is used for the fit. 

11. Comparison of the model results for the Lyapunov exponent (1st iterate) for the entire 
map with the data. The model values, given by a dotted line (simple model) and a solid 
line (improved model), are again indistinguishable. The data point at log 10 (0.128 x 
1(T 3 ) = -3.893 ... is used for the fit. 

12. Comparison of the model results for the Lyapunov exponent for the s = 1.0 simulation 
with the averaged gate data from Fig. 8. The model values are given by a dotted line 
(simple model) and a solid line (improved model). 
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